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ABSTRACT 

The  parametric  decay  of  an  Intense,  coherent  electromagnetic 
wave  Into  a plasma  wave  and  scattered  electromagnetic  waves  In  a 
homogeneous  plasma  (stimulated  Raman  scattering)  Is  generalized 
to  Include  the  coupling  of  secondary  waves  In  the  plasma.  The 
secondary  waves  consist  of  a plasma  wave  harmonic  and  Its  asso- 
ciated mixed  electromagnetic-electrostatic  sidebands.  Explicit 
results  are  obtained  for  wave  amplitudes  and  nonlinear  growth  rates 
when  Yt  ^ 1 , where  y Is  the  linear  growth  rate  of  the  funda- 
mental. It  Is  shown  that.  In  this  regime,  the  nonlinear  growth 
rate  Is  2y/3  . 
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Fig*  1 Trajectories  In  the  x-y  phase  plane.  Circle  corresponds 

to  Initial  values,  x*  .56,  y-4.2,  (J»  ■Tr/3  . Trajectories 
are  bounded  by  curve  x*y  ■ 1. 

Fig.  2 Nonlinear  evolution  of  x and  y . 


Fig.  3 Nonlinear  evolution  of  |P^|  and  |P2kl  arbltary  units 


Fig.  4 Trajectories  In  phase  plane  of  a^*  and  a^^^  . The  curves 
labeled  (+) , (-)  denote  trajectories  for  which  (|)  ■ 0 
and  IT  respectively. 
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I.  INTRODUCTION 

As  is  well  known,  there  are  many  Important  fundamental  processes 
associated  with  the  Interaction  of  an  Intense  electromagnetic  wave  with 
a fully  Ionized  plasma- These  processes  have  Important  Implications 
for  astrophysics  as  well  as  for  laser  studies.  Extensive  investigations 
have  been  made  Into  stimulated  Raman  scattering,  one  Important  process 
In  which  an  Intense  electromagnetic  wave  Impinges  on  a plasma  and  scat- 
ters  off  an  electron  plasma  mode.  If  the  Incident  wave  Is  Intense  enough. 
It  modifies  the  linear  properties  of  the  medium,  thus  creating  many  In- 
teresting new  phenomena.^  In  this  case,  the  Bohm-Gross  mode,  which  nor- 
mally has  a frequency  approximately  equal  to  the  plasma  frequency,  can 
exhibit  a large  amount  of  dispersion.  In  fact.  If  the  Incident  field  is 
large  enough,  we  show  that  electrostatic  modes  exist  with  frequencies 
equal  to  half  the  plasma  frequency,  and  that  they  couple  strongly  to  a 
free  plasma  wave;  l.e. , a wave  which  satisfies  the  dispersion  relation 
describing  a system  In  which  there  Is  no  Incident  field.  In  this  case, 
although  the  back-scattered  waves  associated  with  the  fundamental  are 
several  orders  of  magnitude  larger  than  the  fundamental  Itself,  those 
associated  with  the  harmonic  are  negligible.  Since  processes  which  In- 
volve coupling  of  energy  Into  electrostatic  waves  with  little  back-scatter 
may  be  Important  for  plasma  heating,  we  study  this  process  In  detail. 

We  find  that  this  coupling  results  In  a strong  nonlinear  modification  of 
the  usual  linear  growth  rate  associated  with  Raman  scattering. 
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11.  GENERAL  FORMALISM 

We  assume  that  the  unperturbed  state  consists  of  a large  amplitude, 
plane  polarized,  coherent  electromagnetic  wave 


E = 2E  ^ co8(ic  ‘x-to  t) 
o o o o o 


propagating  In  a homogeneous  plasma,  and  that  the  pair  (a)^,k^)  satisfies 
the  usual  dispersion  relation 


u)  ^ + k ^c^ 


In  the  absence  of  damping.  This  state  Is  characterized  by  electrom. 
oscillating  rapidly  In  the  field  of  the  Incident  wave  E^  with  nonrela- 
tivlstlc  quiver  velocity  |V^|  eE^/mii)^  , while  the  Ions  form  a stationary 

A 

background  . 

Next,  we  perturb  the  equlllbrliim  and  consider  the  mode  coupling 
processes  associated  with  the  medium  In  the  presence  of  the  field  . 

•In  particular,  we  consider  scattering  processes  off  electron  modes  and 


find  that.  In  the  regime  (V^/c)^  - w/uj^  « 1 , modes  with  frequency  of 
the  order  of  half  the  plasma  frequency  exist  and  couple  strongly  to  the 
Bohm-Gross  mode. 

We  examine  the  effects  of  a density  perturbation  associated  with 
an  electrostatic  wave  (u),k)  and  Its  first  harmonic  (2co,2k)  and  look 
for  conditions  under  which  the  harmonic  can  grow  significantly.  Associated 


In  our  formalism,  corresponds  to  a complex  velocity  amplitude. 
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with  these  waves  are  their  satellites  (w  ± jw  i jk.  ) and 

o o 

(2U)  ± joj  ,2k  ± jk  ) where  j is  an  integer.  VJaves  associated  with  har- 
o o 

monies  of  the  pump  , however,  give  contributions  of  order  V^/c 

to  the  current  densities  and  are  therefore  neglected.  We  assume  that 
a)/k  is  much  larger  than  the  electron  or  ion  thermal  speeds  and  therefore 
we  neglect  any  strong  wave  particle  interactions.  Under  these  conditions, 
we  can  use  an  essentially  fluid  model.  All  currents  are  first  calculated 
from  a zero  temperature  fluid  model.  Linear  contributions  to  the  currents 
are  then  modified  to  include  thermal  effects  by  introducting  the  approp- 
riate linear  dielectric  constants  and  susceptibilities.  The  appropriate 
fluid  equations  are  then  simply 


V + (V«V)V  = - - (E  + - V X B) 
dt  me 


n + div  nV  = 0 


(3) 


where  n and  V refer  to  the  density  and  velocity  of  the  electrons.  Ion 
contributions  to  the  perturbed  current  densities  are  neglected  since  they 
are  negligible  in  this  frequency  regime. 

With  this  framework  in  mind,  the  fluid  equations  are  solved  for 
the  Fourier  amplitudes  of  n and  V and  the  relevant  current  densities 
are  computed.  Next,  time-varying  equations  for  the  wave  packet  amplitudes 
are  found,  assuming  that  all  waves  are  coherent.  In  a Fourier  representa- 
tion, Eqs.  (3)  reduce  to 


V - - - Is.  f 

k mO)  mu)  J 

^1,  ~ • 1 n,  V.  dX 

k 0)  J Tti  kj 


0)2 


mo). 


(5) 


3 


where  dX  = dfeidfe26(fe-fei-fe2)  and  dfe^  = dk^do)  (1*1,2).  ’ Labeling 
the  sideband  modes  and  (2w±u)^  ,2^±^^)  by  1,2±  respec- 

tively, and  examining  Eqs.(5)  In  a regime  In  which  ponderomotlve  forces 
are  Important,  we  find  that 


k 'V,  k ‘V 

k . , ^ o 


[k  V n 

— (V  +V  )*E  - (k  E + k E ) • -2. 

(0  + -0  --  ++a)J 

r\  A 


[ (V  + — E ) - (V  - — E )1 
L + mto  + - m<o  - J 


n^  - ± — • (n  ± n,  V ] 
± u)  o ± k o 


V - + — E + li_  r±  )L  (V  ) + 
± mo)  ± mtoLOJ  OK 


— (V,  ‘E  ) 
k o J 


1 -►  le  -► 

+ — k*V  (V,  + — E,  ) 
0)  o k mu  K 
o 


k*V  k*V 

n.,  * n + (n_  -n„^)  2. 

2k  o u '2-  2+  (0 


k Vi 

'2k  - ifc  "2k  - [sr  <^2-^^2+>*"o  - <‘^2-"2--^  V"2-H>  * ifJ 


^ [<v  * S-  - "i--  ^ '2- 

*•0  o 
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^2±  = - izr  * ^V2±  - ^2k 


u - is  t - r 

^2±  “ inu)„  ^2±  ■'■  ina)„  L 


-[ii^  ‘V'o>  * S 


+ + le_  I \ 

^ 0)  r2k  ^ 2inw  ®2k/ 


where  we  have  denoted  the  complex  amplitudes  E.  by  ^ and 

^o  *-o  o o 

neglected  the  effects  of  pump  depletion.  We  have  also  assumed  that 

V , may  be  of  the  same  order  as  V,  . These  equations  can  be  readily 
2k  k 

solved  up  to  terms  of  order  Woi)„  . The  results  are  as  follows: 


n,  * n„  

k o 0) 


V,  = - — P, 
k mo)  k 


n , = n, 

± 0)  k 

o 


t ^ 1 

-►  _ + -►  -► 

V,  = + — [E^  ± — (V  -P,  )] 
± ^ mo)  ± wok 


n-,  “ n 

2k  o w 


V = - P 
’ 2k  2mw  2k 


2k ‘V  K 

"2±  ° "w  ”2k  ’ ^2±  “ + ^ ^^2±  " lia  ^'^o'^2k^^ 

o o 


where 


P^  ”■  2 - — (E  + E ) • V (8) 

k k + - o 

o 

with  a corresponding  expression  for  P2j^  . According  to  Eqs.(7),  the 
ponderomotlve  forces  modify  the  low  frequency  velocity  amplitudes  by 
changing  Ej^  to  and  £211  ^2k  ’ affects  all 

density  contributions.  We  note  that,  in  this  regime, 

of  the  same  order  as  E+  and  cannot  be  neglected.  Therefore,  there  are 
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also  ponderomotlve  modifications  to  the  velocity  amplitudes  associated 
with  the  high  frequency  side  bands.  The  second  order  contributions,  ob- 
tained from  Eqs.(5),  are  as  follows: 


n^2)  . i f(2)  ^ („(2).  „(2)-  1 ^ 

k (i)Lok  - + oJ  k 

”t"’  - ^ • Mi'’  * “i 


^(2)  ^ ^ ^ (k-V  + v' 

± mW^  UJ^  \''k  o / - 03^  o^\  ^ 


I • Mik’  (“i-’  - "j 


7(2)  . .is-  Ik.  (?«>  ^5(2)\ 

2k  2imo  0)^  \'^2+  ^ ^2-  f 


1 /-*-(2)  -►(2)\  •*■■*■  ->1 

k.v^  + v^, 


(2) 

'2t 

.bi.\ 

“2±  L 

n<2)^ 

2k  ''o 

:(2) 

le 

'lo  1^(2) 

2k*V 

O ^(2) 

2± 

mu). 

2± 

“o  ' 2k 

o I 

0)  '^2k 

o 

where 
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k r 

0)  * [“2 


-►1 


k + "2+’'-k-k*”2-’'-ktk  * "-k’'2k 

0 0 

^ "-k-k  V 

0 

* "-k+k  ^2-1 

0 J 

•E  Wl^2i 

-k  ‘2k) 

(’-k-k;*2j-^  (' 

-► 

-► 

2k*®-k)  ( 

V'-k-k  h:-  V" 

-*• 

2k4.k)“2  - ( 

'-►-*■  \ *2+  -►  -► 

k ‘V  1 — ^ V 'P 

,2+  -k-k  ) 2 0 2k 

0 

AktklVl 

0 

!v'2k)-(Hk)“‘ 

(vy  Mv\ )- ("-•'2-)  Uv\)] 


le 

nbu 


(10a) 


r 

N.  “s:*  I 


“2k'^-k±k  ■*■  "at'^-k  ■'■  “ 

o 


-At] 


t 

V± 


*l(*2tAk)+^(jjlAk*»  )] 


♦^[(2^-’-kikJ«2  - (^2.Ak)-r(?,Ak) 
- (H.  )“2  ♦ ( AAk)*t(^,-?k)]i^ 


(10b) 


t 
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“ak  ■ K -[Vk  ”-'J 

'’ik  - - ss  U ''k'^k'  ♦ 5;  «--V  * ir  <»+-t)J 

+ ^ [(S-\)kp  - (k.-v^)E.(v^-Pt)  - <2^,-f.)J+(v„-p^)]  (loc) 

"2t  - ^ • Mt  * "t\] 


f.  is-Fi^ 

'21 


+ --^[-(^  ‘V,  )k  (V  •?  ) + (it'V  )fcl 

(1)2+  L ± k'  ±'  o k'  ± J moxjjjj 

P - (E  +E  ) -V  , P - (1  .+1  ) 

+ “ O 2 2+2-0 


(lOd) 


We  now  apply  the  theory  to  a Stokes  mode,  D «D^  , where' 


- k^^c*  - 


e.  - 1 - 0)  ^/i)ij 
± P ± 


This  implies  that 


"+  ~ \ ~ f ^- 


(13a) 


ht  " c ®2k 


(13b) 
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as  we  show  below.  With  this  ordering,  Che  second  order  perturbations 
reduce  to 

/k*V  \ / k»V  \ . k. 


„(2)  (2) 

"t  i— K 


V^2)  . ^ /y  ,y  \ 

k 0)  ^'^2k  '^-k'' 

(14a) 

^(2)  _ ) 
\ 0)^  /'  2k  -k^ 

0) 

(14b) 

..2 

V^2)  , ^ _k_ 

'^2k  2o) 

(14c) 

2 ^ 

V<2)  . k/  oW 

2±  2o)\  0)  /‘"2± 

(14d) 

The  current  density  perturbations,  due  to  the  response  of  the  electrons 
v'.o  the  electric  and  magnetic  fields,  are  then  as  follows: 


2 

0) 

_E 

-► 

P . 

e 

. 3 

2 

0) 

-E.  k 

0) 

4tto.’ 

2 

0)2  m 

2 

0) 

JE 

\K 

± ic  1 

V ‘P  V 

L ± 

±\ 

) w 

/ o 

(15a) 


0)  kP, 


i — — -E.  ^ p p 

4Tr  2 u*  m u)^  2k  -k 


2 

P.  0) 


► 4 n-^  «a 

I * P + II  • Jl  V -ii  — t. 

2k  4TT  20)  2k  2 m 


(15b) 


(15c) 


2±  4tt  I 2±  2u)  o *^2k  o)  *^2k  2o)  oj 

^ - 0)*  / kP.  \2 

± — -I  V — ^ 

•iTT  2 0)2  o V mo)2  / 


(15d) 
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Thermal  corrections  to  the  first  order  current  density  perturbations  can 

2 9 2 

be  obtained  by  letting  -*•  -X(lt,u)  , u /(D . etc.  The 

P P - - - 

second  term  on  the  right  side  of  Eqns  (15b)  and  (15d)  will  be  shown  below 
to  give  a negligible  contribution  to  Eqns  (17a)  and  (17b),  respectively. 

If  we  add  thermal  corrections  to  all  other  terms  of  first  order,  we  obtain 
the  following  expressions  for  the  first  order  current  density  perturbations: 

t(l) 


(16a) 


r(i) 


ICJ 


7T  ‘^+X+E.  ? X(k,u))V  VV.  t r-E-  kj 

ATT  ±''±  + ATT  ’ o k ATIO)^  ±\  0)  / 


■ - S 2“X(2k.2a.)5j, 


(16b) 


(16c) 


►(1) 

’2± 


Att  ‘^2±^2±*^2±  ^ ATI  X(2k,2a))V^2kP 


lh)p 


_ir  (V!2k) 

'2+  2±\  20)  / 


o 2k  ■ Attu), 


(16d) 


Since  the  second  order  current  density  perturbations  ^2k^ 

both  proportional  to  k , the  fields  and  remain  electrostatic 

In  second  order. 

The  Fourier-transformed  wave  equations  for  E^  and  £2^  are  as 


follows: 


2 2 -K  ^ 

9 0 ^ ^ ^ -t.  W /V  ‘P,  \ 

|(k.-(.,/c»)I  -kjkJ.E,  -Jf  XtE,*  JT  xVk  * 


2 

~ 0).  0)  V„  , 2 

^ ~ le  — ^ — E-  — — p p 

2 c2  0)2  m 0)*  *^2k*^-k 


(17a) 


I 
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1 


I 


[(k2i  - - '^2±^2±’*^2±  ' h±h±  - ^ X(2k,2u,)V^2kP2^ 

u»^  ’PouX 

^ c2  ‘‘2±\  2ui  f 


- V.2* 

+ 1 ie  2±  _£  _o  k_  2 

2 c*  U)2  m (1)^  K 


(17b) 


where  I is  the  unit  operator.  These  equations  reduce  to 


[D^I  - k^k^c^l'E^  =*  V^Pj^  + 


^*^2±i  ~ ^2±^2±^  ^*^2±  * '^2±^2k  ^ ‘^2± 


where 


V - V ^ (D  k. 

1 p ± ko) 


(20a) 


V = 0)  xV  k 
o^  o 


(20b) 


and  with  a similar  expression  for  V2^  • The  normalized  current  densities. 


, and  are  given  by 


= ■ I i#  ^ Pzi^P.k 


(21a) 


^2±  = * I ^ Pk 


(21b) 


where  terms  of  order  w/wq  are  neglected.  Similarly,  Ej^  and  Ezk 
satisfy  the  eqiiations 
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(22a) 


(2u>)*E2k  = - (2u))^x(2k,2a))P2k  + 


(22b) 


where 


4k  * ■ 5^  ^ ’’k  • 


We  now  invert  Eq,(L8)  and  obtain 


- ^±*<Vk  + 1+) 


where 


II  - k+^+cVu)*e^]/D^ 


(23a) 


(23b) 


(24a) 


(24b) 


Using  Eq8.(22)  and  (24)  to  eliminate  E^  , E^  and  E_  from  Eq.(8),  we 


obtain  the  following  equation  for  : 


**lr  ^ X('^.w)  + 


K * - 8-  *-))] 


4 _ .M.  / 

(0  ^ 


(25a) 


(25b) 
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We  note  that 


Vq-  ’V 


U)  V • (>r - M -k  ) 

p ku  o -f  + — 


and  therefore,  that  Eq.(25)  reduces  in  lowest  order  to 


Ji  kV 

P A((jj,k)  - ^ 

k ’ 0)2  (4) 


• M • 


where 


A(o),k)  - 1 + x(k,o))  + k^x(k,o))V  .(M  + ) -V  (28) 

O "T"  — o 


Is  the  modified  dielectric  consteint  of  the  medium  in  the  presence  of 
■*  * 

£_•  Referring  to  Eqs.(21)  and  (23),  we  see  that 


kV  0) 


^V.  O'  0) 


which  allows  to  be  eliminated  from  Eq.(27)  to  give 

Pj^A(o),k)  ’ ^ [l  + k^V^-(tf^  + MJ.vJ 


0)*x(k,0)) 


(-1  + A(o).k)) 


Eq. (28)  can  be  readily  shown  to  reduce  to  the  result  of  Drake  et  al.. 


if  one  remembers  that  V ■ - ief  /mo)  , while 

o o o 


in  Ref. 5,  V = — . 

o mo) 

o 


t 
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where  £q.(28)  has  been  used  to  further  simplify  the  result.  In  a similar 
manner,  we  find  that  the  equations  for  and  A(2a},2k)  are  the  same 

as  Eqs.(30)  and  (28)  respectively  if  (u),k)  is  replaced  by  (2u),2k). 
Using  Eq. (30)  to  eliminate  from  Eqs. (22)  and  (2A,  gives 


'A(a),k)Ej^  - 


a)^A(u>,k)E^  ' - 


Similarly,  we  find 


(2a))2A(2a),2k)E2^  • 


(2a))^A(2a),2k)E2^  - - 2(A)^k  ^ (34) 

2 2 

For  a Stokes  mode,  D ^ ukjO  , ~ oj  and  O)  for 

- 0 + 0 2±  o 

k = 2k^cos6  .*  Estimating  the  relative  magnitudes  of  Ej^  and  E^  from 

Eqs. (31-32),  assuming  that  (V  /c)*  ~ w/o)  , yields 

o o 


°+  , gjo  1 

^ '"o'" 


E u k 1}  'V  u)  kV 
0-0  o o 


(35a) 


(35b) 


Similarly,  we  find 


0)  kV 
o o 


We  now  show  that,  for  ui  - , A(2a),2k)  ■ 0.  According  to  Eq.(28), 

A(2aj,2k)  - e(2a),2k)  + (2k)2x(2k,2a))V^*(tJ2+  + M2  > (37) 

since  e ■ 1 + X • The  last  term  on  the  right  side  of  Eq.(37)  Is  of 
order  (V^/c)^  , while  G((i)^,2k)  = 0 . Therefore,  A(2t*),2k)  - 0 for 
some  0)  close  to  . 

If  A((i)./2,k)  « 0 for  some  E , there  may  be  significant  coupling 
X o 

of  energy  between  the  waves  and  E2j^,E2^  . The  dispersion  rela- 

tion A(ciJ,k)  “ 0 for  a Stokes  mode  In  the  regime  (i)»kV^  reduces  to® 


(oj-Aw)  (0)^-0)^^  +1(1)^  ImX) 


- 2a)o‘^p''^'(|Vor/c^) 


where  the  geometry  of  the  scattering  is  described  by  the  following  para- 


meters: 


slni|)  COS0 


sln^cj)  “ |k  X V I ^/k^  |V  | ' 
' - o'  - ' o 


Aw  “ c 


k*k 

2 o 


D (k  ,0)  ) = 2u}  (w-Aw) 
— — — o 

k = 2k  C089 
o 


We  neglect  Imx(k,u))  compared  to  Imw/Rew  and  look  for  unstable  modes 
with  Re  0)  = . Let 


F(a))  ■ (u-Aw)  (w*-u)|j^) 


We  plot  F(b))  vs.  b)  for  real  values  of  u and  determine  Aw  such 


i 

fc 
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that  F'(u)^/2)  - 0 . This  requires  Aw  * - ^ . For  this  value  of 

Aw  , F(Wj^/2)  ■ ~ • Eq.(38)  has  three  real  roots  for 

A ■ 2w  w ^ l*/c^)  < W-’  . For  larger  values  of  A , two  real 

op  o lo  X. 

roots  disappear,  and  the  system  becomes  unstable,  there  being  now  one 
real  root  and  two  roots  which  are  complex  conjugates  of  each  other. 
For  back-scattering  , = 1 , which  Implies  that 

^32  IT 

o 

Is  necessary  for  instability  at  this  frequency.  Setting 


w - w^/2  + lY 

A • W.^ 

o 16  £ 


Aw  ■ - T-  w.  + Aw 

H £ 


and  assuming  that  Imy  “ 0 , the  real  and  Imaginary  parts  of  Eq.(38)  take 
the  form 


Wj^  - Aw 


y - ! “( 


* AA 


(A3a) 


Aw  w^  + y*  - 0 


(43b) 


where  AA  ■ A - A^  and  the  relevant  solution  for  Aw  , 


4w  - f (-^[l  -Vl  + H v] 


can  be  expanded  to  give  the  following  results  for  Aw  and 


i 

• 0* 


Aaj 


2 AA 
5 (HP 


- 2 M 

In  terms  of  the  pump,  we  have 


(45a) 

(45b) 


Q AE 

o 

where  is  the  threshold  value  corresponding  to  Y ■ 0 and  AE^  is 

a small  Increment. 

Applying  standard  techniques  to  Eqs. (31-34)  enables  us  to  obtain 
the  following  results  for  the  equations  describing  the  wave  packet 
amplitudes,  assuming  that  they  are  spatially  independent:^ 


3(ot)^A).  ‘^k 


3<i) 


k/k 


-w  kM^*V 
o ~±  o 


(47a) 


(^*\) 


"2k 


"2± 


(w^A) 


3oo 


k/k 


-40  2kM,.*V 
o ^2±  c 


(47b) 


where 


Afti  ■ u)(2k)  + u)(-k)  - a)(k) 
Aflz  - 2a)(k)  - aj(2k) 


i /-LV  J_ 

“ 3 la'  J “l  “ 15  E 

Xj  O 


(48a) 

(48b) 

(48c) 


u«0j^/2 
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(48d) 


Since  A(2u,2k)  - c(2u),2k)  , Yt  just  the  usual  linear  Landau  damp- 

1j 

Ing  rate  of  a Bohm-Gross  mode  and.  In  this  case.  Its  sidebands. 

If  we  neglect  the  slight  frequency  mismatches,  and  An2  > 

the  equations  for  the  wave  amplitudes  reduce  to 


■ - <=■“<,"  «t-’o  V-V 

idi  * *’k 

(w  \)*2i‘'>  ■ «“o"‘  »2t-’„  ■’k' 


where 


Cl 


1 ^ /x  y 

2 mo)^  / 


C2 


, ek 
mw^ 


(49a) 

(49b) 

(49c) 

(49d) 

(50a) 

(50b) 


S 


\ 

\ 

\ 
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III.  ANALYSIS 

We  now  discuss  the  properties  of  Eqs.(49).  We  assume  that  at 
t*0,  ~ Ej^  . Then,  according  to  Eqs.  (35-36),  E^  ~ (V^/c)Ej^  , 

E_  ~ Ej^/(V^/c)  , and  £2^  (V^/c)Ej^  . Since  a)»kV^  , ' 

We  show  below  that  E2j^  and  E2^  vary  on  a time  scale  ~ l/y  . 

We  therefore  neglect  Y^  Eqs.  (49).  Also,  a knowledge  of  , 

E2j^(t)  and  the  Initial  conditions  of  all  the  waves  Is  sufficient  to 
determine  the  nonlinear  evaluation  of  E^  and  E2^  . One  merely  ex- 
presses Eqs.  (49)  In  Integral  form  as  follows: 


E^(t) 

E^(t) 


ft 


dT  e^^‘''‘'^,P2^(T)P_j^(T) 


E.  (0)e^*^  - 0)  k(M. ‘V  ) 
± o ~±  o 


dT  e^^‘'"’'^CiP2i^(T)P_^(T) 


(51a) 

(51b) 


with  similar  expressions  for  E2j^  and  E2^  . Solving  for  the  Integral 
In  Eq.(51a)  allows  E^(t)  to  take  the  form 

E^(t)  - E^(0)e^*'  - ^ok(Mi*Vj,)(Ej^(t)  - ^(0)6^*")  (52) 

The  equation  for  E2^(t)  Is  obviously  the  same  with  Y“0  and  ((0,k) 
replaced  by  (20), 2k),  Instead  of  analyzing  ^(t)  and  Its  harmonic.  It 
Is  more  advantageous  to  examine 

Fourier  amplitudes  of  the  total  first  order  field  that  govern  the  motion 
of  the  electrons.  In  terms  of  P^^  and  E2k’  find  that 


c,P2tP-k(-l  + A(a),k)| 

0)=0)j^/2 


(.53) 
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and,  similarly. 


at  ^2k  “ *^2^^ 


where  we  have  taken  appropriate  linear  combinations  of  Eqs.(49)  and  used 
Eqs.(28)  and  (37),  remembering  that  A((Dj^/2,k)  = A(a)^,2k)  ==  0 . We  have 
therefore  reduced  the  analysis  of  a six-wave  system  to  an  analogous  one 
Involving  only  two  waves. 

In  order  to  analyze  Eqs. (53-54),  we  first  transform  them  to  a dim- 
ensionless form  and  then  express  them  In  a polar  representation.  In 
terms  of  the  nondlmenslonal  quantities  r * yt  » ^ “ (/cjCj/y)?]^  » 

A2]^  = ’ Eqs.  (53-54)  take  the  form 


( 9t  “ “ " ^2k\ 


(55a) 


— A = A.  ^ 
3t  *2k  \ 


(55b) 


where,  denoting  the  complex  conjugate  by  * , , since 

Is  a reality  condition  on  the  total  electric  field.  Setting 
Aj^  ■ a^e^*^^  and  A2j^  “ a^j^e^'^^k  ^ Eqs.  (55)  reduce  to 


dx  *k  ^ 


- a2^a^  cos 


(56a) 


^®2k“  V«8<(. 


(56b) 


^ (V®2k  ‘ »2k  ♦ 


(56c) 


where  <t)  - . Eq.(56c)  Is  obtained  from  a linear  combination 

of  the  Imaginary  parts  of  Eqs.(55a)  and  (S5b).  Its  solution 


®k*®2k  “ ®k^  (0)a2j^(0)  aln  ()i(0)e 


allows  Eqs.(56a)  and  (56b)  to  be  further  reduced  to  the  form 


^ + 2/x^y-l  a 


(58a) 


^ - 2/x^y-l  e^^  a 


(58b) 


where 


x(t)  - a,^(T)( 


.-2t  -Vs 


(59a) 


2 yu 


y(T)  - ® 


(59b) 


a ■ a^  (0)a2j^(0)8ln  (j)(0) 


(59c) 


by  allowing  cos  ()>  to  be  eliminated.  We  note  that 


dx  -2t 

dy  “ ■ ® 


a relation  which  enables  Eqs. (58)  to  be  readily  analyzed  In  the  phase 
plane.  An  upper  bound  may  be  obtained  for  x(t)  by  Integrating  Eq.(60). 


We  find  that 


-2t'  -2t' 

dx  - [-d(e  y(T'))  - 2e  y(T')dT'] 

■'x(O)  ^t’-O 
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which  reduces  to 

x(T)  s y(0)  x(0)  (62) 

where  the  Inequality  is  a result  of  neglecting  negative  definite  terms 
from  the  right  side  of  Eq.(61).  We  show  below,  In  fact,  that  x(t)  goes 
as3rmptotlcally  to  zero. 

It  is  useful  to  examine  the  motion  of  x and  y In  the  phase  plane. 
According  to  Eqs.(58},  the  motion  Is  bounded  by  the  curve  x*y-l  ■ 0, 
and  from  Eq.(60),  |(dx/dy)|  Is  equal  to  one  at  t»0,  and  decreases  mono- 

tonlcally  to  zero.  The  motion  can  be  analyzed  further  by  examining  an 
equation  which  describes  the  evolution  of  (|)  . We  take  an  appropriate 
linear  combination  of  the  Imaginary  parts  of  Eqs.(55)  and  get 

y ^ * (-  xe^^+2y)a  ^*sin  (|)  . (63) 

Turning  points  in  Fig.  I occur  when  (p  ■ 7r/2  or  37r/2  . Since 
|dx|  s |dy|  , between  turning  points,  the  trajectory  to  the  right  will 
always  lie  above  the  last  trajectory  to  the  left,  and  also  above  the 
trajectory  to  the  left  that  Immediately  follows*  The  motion  Is  relatively 

ii 

rapid,  depending  on  the  magnitude  of  a (nominally  In  the  range  1-10), 
except  near  the  boundaries.  As  a first  approximation,  we  assume  straight- 
line  motion  between  turning  points  except  near  the  boundaries  (see  Fig.  1). 

We  define  t%K>  sets  {x^  } and  ^ such  that  labels  the 

value  of  X at  a left  turning  point  and  x^^  , one  that  occurs  to  the 
right,  where  It  is  understood  that  j Increases  with  time.  It  is  clear 
from  the  analysis  above  that  both  sets  are  monotonlcally  decreasing  with 
lower  bound  equal  to  zero.  Therefore,  x(t)  goes  to  zero  asymptotically. 
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Nonlinear  growth  rates  may  be  determined  from  an  analysis  of  Eq.(63). 
First,  we  note  that  the  signed  square  root  In  Eq8.(58)  Is  proportional 
to  cot  <|>  . According  to  Eq.(S7),  If  (<>(0)  j*  0 or  it  , then  4i(t) 
never  goes  to  0 or  tt  In  finite  time,  unless  perhaps  x or  y 
becomes  Infinite  at  that  time.  The  motion  Is  therefore  restricted  to 
values  In  the  half  circle,  0<(^<ir  or  7T<4)<27r.  We  assume  for 
definiteness  that  0 < 4)(0)  < ir  , and  note  that  Is  equal  to 

zero  somewhere  an  every  trajectory  between  turning  points.  Let 
{Xj-  x(Tj)}  be  the  set  of  values  of  x such  that  (|)'(Tj)  « 0 . 
According  to  Eq.(63),  at  this  value  of 

- x(Tj)e  ^ + 2y(Tj)  - 0 . (64) 

Not  only  does  x(Tj)  go  to  zero  as  J , but  also  dx/dx  goes  to 
zero  with  Increasing  j as  exhibited  In  Fig.  I.  This  Implies  that 

x^(Tj)y(Tj)  - 1 -<■  0 (65) 

asymptotically,  as  we  see  from  Eqs.(58).  The  envelopes  of  the  curves 

describing  x(t)  and  y(T)  are  now  readily  determined  by  examining 

their  values  at  the  turning  points.  After  several  trajectories  between 

turning  points,  the  motion  will  be  nearly  horizontal,  so  that  slight 

2t 

changes  In  the  slope  of  the  line  -xe  +2y  - 0 will  Intersect  the  curve 
x^y-1  at  neighboring  turning  points.  Therefore, 

x*y»l  , x ■ 2e  ^^(l+6)y  , 0 < 6 « 1 (66a) 

(at  a left  turning  point) 

while 

x^y-1  , x - 2e"^'^(l-6*)y  , 0 < 6'«  1 (66b) 

(at  a right  turning  point) 


t 
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from  which  we  find 


Max  X ■ 


n r 

Max  y - y.  - Tr 

TSTj’^  J [4(1-6’)1  " 


or  in  terms  of  |Pj^(  and  1^211'  ’ 


Max  I\1  * 7=  a'^[2(l+6)l 

T St j ^ /ciC2 


^ Vs  Vst/ 

Max  |P,.  1 - ^ ^ V 

TsT  r 2k.  Cl  [4(i.fi.)2]76 


(67a) 


(67b) 


(68a) 


(68b) 


Max  iPj 
Max  fpT 


■V  ^ -VI  iv  ^ 


X , r , . 


Tj"  > Tj  (68c) 


Since  T = Yt,  this  implies  that  Pj^  and  P^j^  evolve  nonlinear ly  with 

2 

a growth  rate  Yjj^  " 3 • 

The  case  where  <^(0)  - 0 or  sr  must  be  considered  separately. 
According  to  Eq.(57),  this  Implies  that  4>(t)  ~ 0(0)  for  all  T . In 
that  case,  Eqs.(56a)  and  (56b)  reduce  to 


F *2k  - * *k' 


(69a) 


(69b) 


(69c) 
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from  which  another  constant  of  the  motion  may  be  deduced,  namely. 


% 


aj^^T)  + a2^(T)  + 2a2j^(T)  = (70) 

2 

A typical  trajectory  (aj^  vs.  a2k^  shown  In  Fig.  4.  The  orbit 
labeled  (+)  denotes  ())(t)  = 0 , while  the  orbit  denoted  by  (-) 
refers  to  (|)(t)  = it  . The  small  circle  corresponds  to  the  Initial 
point  , T = 0 . According  to  the  figure,  both  orbits  touch  the 
a2jj  axis  and  stop  at  a time  which  can  be  shown  to  correspond  to  a 
few  units  of  T . In  this  case,  therefore,  Pj^  0 , and  P2j^ 
saturates  after  a few  linear  e-foldlng  times  (yt  ^ 1)  . (Negative 
values  of  a2j^  on  the  orbit  labeled  (-)  are  due  to  the  fact  that 
the  trajectory  of  A^^  In  the  complex  plane  Is  along  a ray,  starting 
at  a point  In  Quadrant  I , passing  through  the  origin,  and  ending  In 
the  third  quadrant.) 


IV . CONCLUDING  REMARKS 

We  have  examined  the  nonlinear  properties  associated  with  stimu- 
lated Raman  wave-wave  scattering  off  a natural  plasma  mode  and  a 
strongly  modified  plasma  mode,  and  found  that  this  complicated  process. 
Involving  the  mode-coupling  of  seven  waves,  can  be  reduced  to  an  analy- 
sis of  a two-%»ave  system.  An  Interesting  feature  of  this  system  Is  that 
the  waves  have  no  asymptotic  dependence  on  initial  conditions,  other 
than  an  assumption  that  E|^  and  its  harmonic  are  initially  of  the  same 
order.  We  can  then  determine  analytically  the  asymptotic  prop- 
erties of  the  waves,  and  find  that  there  Is  a strong  nonlinear  modi- 
fication of  the  usual  linear  growth  rate  associated  with  stimulated  Raman 
scattering  and  also  that,  although  the  harmonic  has  a rather  small  linear 
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damping  rate  y , it  grows  nonllnearly  at  the  same  rate  as  Its  funds- 
mental  in  the  asymptotic  limit. 
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